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Malmquist bias is present in all astronomical surveys where sources are observed above an 
apparent brightness threshold (e.g., Malmquist 1920; Tammann 1987). Those sources which can 
be detected at progressively larger distances are progressively more limited to the intrinsically 
luminous portion of the true distribution. This bias does not distort any of the measurements, but 
distorts the sample composition. 

We have developed the first treatment to correct for Malmquist bias in linear regressions of 
astronomical data. This poster presents a demonstration of the corrected linear regression that is 
computed in four step and illustrated by the two figures below. 

Step 1: Create a synthetic parent sample of galaxy luminosities in the blue (Lb) and 40-120 pm 
infrared ( Lfir ) for which the intrinsic luminosity dispersion and correlation are known. 

Step 2: Place the parent galaxies at randomly selected distances and observe them with a detection 
threshold, thereby creating a data sample with Malmquist bias. 

Step 3: Compute the linear regression of the data sample both with and without correction for 
Malmquist bias. 

Step 4: Performs steps 2 and 3 repeatedly, generating two distributions of regression results. It 
will be seen that the peak of the corrected regression fits is much closer to the true functional 
relationship of the parent population. 

In Verter (1992) the corrected regressions are applied to the same galaxy sample that was 
previously studied by both Lonsdale-Persson and Helou (1987; PH) and by Trinchieri, Fabbiano, 
and Bandiera (1989; TFB). We show that once Malmquist bias is removed, all the luminosity- 
luminosity plots agree, to within 2 a , with a pure increase of galaxy luminosity with galaxy size. 

Figure 1 is a parent sample of 200 Lb and Lfir values generated by Monte Carlo simulation. 
It is assumed that log(La) is described by the Gaussian distribution fit to the Virgo cluster by 
Sandage, Binggeli, and Tammann (1985), and that \og(LpiR/ L b) is described by the Gaussian 
distribution fit to UGC galaxies by Bothun, Lonsdale, and Rice (1989). The parent population has 
a built-in linear correlation between Lb and Lfir- 

Figure 2 is one of the possible data samples generated by placing the parent galaxies at var- 
ious distances and observing them with a surface brightness detection threshold. The observing 
parameters were chosen to match the previous study by TFB; the galaxy distances are uniformly 
random between 12 - 60 Mpc, and the detection thresholds are Bj < 13.2 mag and / fir > 1.5 Jy. 

The solid fine in Figure 1 is the best fit of an ordinary linear regression; it has slope 1.028. 
The parent population was simulated only once, but the observed data sample was simulated 450 
times; the number of detected galaxies was typically between 35 and 55. The dotted line in Figure 
2 is the slope of the ordinary linear regression fit to the data samples, 1.386iouo3- The dashed line 
in Figure 2 is the slope of the weighted linear regression fit to the data samples, 1.1461 q; 296- The 
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lower error bar on the ordinary regressions is too high, by 3 ct, for the fit to be consistent with the 
linear correlation of the parent sample, whereas the weighted regressions are consistent within 1 a. 

Our method of correcting for Malmquist bias uses a computer program that computes the 
linear regression with errors in both variables, but instead of the error in each observation we sub- 
stitute the volume sampled in its measurement. This volume ( V \ ) is determined by the distance at 
which the galaxy surface brightness would fall to the detection threshold. Treating the sampling 
volume as if it were an error implies that galaxies at larger distances receive less weight in the 
linear regression, because they are more subject to Malmquist bias. This is illustrated by the sizes 
of the symbols in Figure 2, which are proportional to log(l/VgVp , /j^). Figures 1 and 2 intuitively 
illustrate why the weighted regression is an improvement: Malmquist bias introduces a dispropor- 
tionate number of high luminosity pairs into a galaxy sample. An ordinary linear regression is 
pulled up towards higher slope, but when the observations are weighted the fit favors the more 
completely sampled galaxies. 
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